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I. INTRODUCTION 


A. BRIEF BACKGROUND AND PROBLEM STATEMENT 

The development of computer vision is an ongoing area of research in the field of 
artificial intelligence. The ability of the computer to identify objects in images, or in 
the real world, opens up a whole new form of data input. Research in the field varies, 
from vision in controlled environments, such as industrial assembly lines, to 
identification of objects in uncontrolled, real world environments. 

Vision research has been broken down into two phases. The first phase, early 
processing, has seen a great deal of experimentation in using the statistical properties 
of texture to identify regions in images [Ref. 1,2,3]. This research has produced many 
complicated mathematical equations to differentiate between textures [Ref. 1,2,3]. The 
purpose of this study is to identify two textures using relatively simple equations. 

In the second phase of computer vision, actual objects are identified. This 1s 
accomplished by matching descriptions. This study does not go into the second phase 
of vision. 

To achieve the desired results of this thesis, it is necessary to first have some way 
of storing an image in the computer. The author along with others developed a system 
for this purpose prior to this study (Ref. 4]. This svstem digitizes and displays images 
up to 1024 by 768 pixels. 

For data, the U.S. Army provided aerial photographs of Fort Hunter Liggett, Ca. 


The objective of this thesis was to pick out road surfaces on these photographs. 


B. METHODOLOGY 

In beginning this study, it was necessary to make a guess as to What qualities on- 
road and off-road regions possess to enable the human vision system to distinguish 
them. Once the qualities were determined, a way of representing those qualities to the 
computer was required. A mathematical strategy was then needed to represent them in 
a simple way. The next step was to process a series of test images to evaluate the 
equations. Finally, it was necessary to analyze the results of the test to determine if the 


equations worked, and which worked the best. 


C. ORGANIZATION 

A brief review of the basics of computer vision and previous work is presented in 
Chapter Two. Chapter Three consists of a detailed problem formulation. Chapter 
Four details the implementation of the recognition system. Chapter Five is an 
evaluation of the results obtained in scanning an image using different methods. The 
final chapter, Chapter Six contains conclusions and recommendations for further work 


in this area. 


II. SURVEY OF PREVIOUS WORK 


A. INTRODUCTION 

In this chapter, a brief background of previous work in computer vision is given. 
The emphasis is on texture analysis, with some description of early vision techniques 
and a brief description of higher level processing. 


B. GENERAL APPROACHES TO COMPUTER VISION 

The general approach to computer vision is to break the process down into two 
distinct steps. In the first step, the goal is to get useful information from the original 
image. The second step tries to identify objects from the information found in the first 
step. Such processing 1s called early processing. In early processing, every part of the 
image 1s processed in parallel. Subsequent processing uses the data produced in early 
processing sequentially and extracts a symbolic description of the image. [Ref. 5: page 
95] 

The way in which computer vision is broken down is similar to human vision. 
Humans appear to make some discriminations in an image automatically and in 
parallel. These early discriminations include texture segregation. In the second stage of 
human vision, focused attention is required and this appears to be performed serially. 
This second stage identifies the objects we see. [Ref. 6: page 156] 

1. Early Vision Image Processing 

In early vision processing, the first step is a description of the edges and other 
interesting features of an image. This description 1s commonly called the primal sketch. 
To produce the primal sketch, it 1s necessary to analyze the image and abstract the 
changes in gray levels to produce a database. This database contains data structures, 
each describing a feature of the image. Each item in the image is made up of one or 
more edgelets, or small edges. Each feature of an image is a combinations of these 
edgelets and are categorized as follows: edges, which consist of edgelets with the same 
orientation, bars, which are two rows of edgelets going in the same direction, and blobs, 
which are edgelets enclosing a small region. [Ref. 5: pages 98-101] 

After developing a primal sketch, there are several ways to extract more 
information. These processes are; motion analysis, which is identifying features that 


seem to be moving together, stereo disparity, which analyses the difference between the 
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two views to determine depth and orientation, photometry, which determines the shape 
of surfaces from shading and grouping processing, that organizes smaller features into 
larger ones [Ref. 5: page 99]. 

Two types of grouping processing that are particularly interesting to the work 
in this study are virtual lines and texture analysis. In the virtual lines technique, the 
computer draws connecting lines between blobs or endpoints. The decision as to which 
items to connect comes from a histogram that gives greater weigh to near objects. In 
texture analysis, the computer tries to identify repeated elements or patterns. There are 
several methods used to accomplish this and they are discussed below [Ref. 5: pages 
Het, 120). 

2. Late Vision Symbolic Description 

The second step in computer vision is scene description. This step is broken 
down into three stages. First, the data produced in the first step, is converted to a 3-D 
object-oriented stage. In the second stage, the computer compares the objects in the 
image with a library of known object types. The last stage is to match an object with a 
known object type, in this way identifying the object in the image. To accomplish stage 
one, a shape notation is required. A shape notation is a standard way of describing 
objects both in the library and in the image. In the third stage algorithms are required 
that take into consideration that an object in an image can be partially hidden, rotated, 
seen from a different angle, or that the object varies from the ‘standard’ object in the 
library [Ref. 5: pages 144-155]. 


C. IMAGE SEGMENTATION 
1. Types of Segmentation 

There are many segmentation methods, with many variations. Three of these 
methods are split, merge, and split and merge. In the split method, a region is 
obtained from the image. This region is compared to known classes. If this region is 
identified as containing only one class, it labeled as that class. If not it is split and the 
process is repeated on each part until all parts are labeled. [Ref. 2: page 274] 

In the merge method, a small area is chosen and classified. Then each 
surrounding area is examined. If an adjoining area is the same class, it is combined 
with the first area and the process is repeated. In the split and merge method, the split 
method is used until a region is classified. Then the merge method is applied 
[Ref. 2: page 274]. 
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Another method is that of using regions called local windows. In this method, 
the image is divided up into many windows. The size and shape of these windows 
depends on the image and the size and shape of the objects to be identified. The size 
and shape of the windows is a critical factor in the process of correctly identifying the 
objects in the image [Ref. 1: page 368]. 

2. Features for Classification 

The regions described above normally are classified in the early stages of 
vision processing using either gray levels or textures. Within the area of texture, there 
are two major analytical approaches, statistical and structural [Ref. 3: page 336]. The 
structural approach is used if the texture has a clear pattern and can be matched with 
stored patterns. The statistical approach is used when pattern matching is not 
possible. 

There are many possible statistical approaches that have been tried, some with 
very good results. In a project at Louisiana State University, six texture statistical 
measurements were chosen [Ref. 2]. These measurements were: inertia, cluster shade, 
cluster prominence, local homogeneity, energy, and entropy. Their results of classifying 
aerial photographs into nine classes produced a result of 90% overall correct 


classification [Ref. 2: page 274]. 


D. IMAGE SEMANTICS AND SYNTACTIC PATTERN RECOGNITION 

In general,the syntactic features of roads are that they consist of two continuous 
edges with a separation not greater then X, and not less than Y. Another feature is 
that these edges or roads are connected; 1.e. a road does not have two ends. Roads can 
have curves, but these curves have a radius no less than Z. Of course the values for X, 


Y, and Z depend on the scale of a particular given image. 


E. HARDWARE AND SOFTWARE SYSTEMS FOR IMAGE ACQUISITION 
AND PROCESSING 


1, Cameras 
a. Pictures 
There are several media for obtaining images. These include slide 
transparencies, cine film, videotape and hard-copy representations such as photographs 
[Ref. 7: page 38]. Most of the work examined in this field uses black and white aerial 


photographs. 
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b. Digitizing 
In order to get an image into the computer, it is scanned, digitized and 
stored in the computer pixel by pixel. Each pixel is identified by an x and y coordinate 
for its location and a gray level if a black and white photograph, or three numbers 
giving the red green and blue intensity of the pixel in the case of a color photograph. 
2. Displays 
Although it is not necessary for the computer in its processing to have a 
display device, it is preferable for humans, because of our limitations, to be able to see 
the results. There are many different graphics devices to choose from. The most 
common display device for interactive graphic 1s the raster scan, which can display a 
mono-color, gray level or color image. 
3. Languages and Computer Systems 
In most research, after the hardware is decided on there is a choice of several 
programming languages. These choices range from assembly languages to high level 
artificial intelligence languages such as Prolog and LISP. Of the languages available, 
Common LISP is particularly appealing because it is a DOD standard. 


F. SUMMARY 

Much of the prior work in computer vision thas been focused on the areas of 
early vision analysis, and in texture analysis. The way vision has been approached is 
parallel to the way it 1s believed that humans do vision processing. The use of texture 
has proven to be very successful in identifying regions of known classes. The 
identification of objects is achieved by creating descriptions of objects and then 
matching those descriptions to known objects. Although research in this field has met 
with success, there still remains much work to be done in both early and late vision 
processing. 
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Ill. DETAILED PROBLEM FORMULATION 


A. INTRODUCTION 

This chapter describes the methods considered in solving the problem of finding 
roads in aerial photographs, and the methods chosen for solving this problem. It 
includes the data format of the files produced by the digitizing camera on the 
VAX/VMS system, and the files produced on the IRIS system containing the scanned 
image road identification. Also described are the features selected for analysis and the 


utilization of these features in identification of on-road areas. 


B. DATA FORMAT 
_ 1. Digitized Image File 
The map data consists of 500 by 500 pixel black and white digitized images. 
This data is contained in a file with a eighty-six byte header followed by the digitized 
image. The header consists of single byte character data and two byte integer data, 
with the least significant byte (LSB) first. The header consists of the following: 
1. The integer 1, which indicates a black and white image. 
2. An integer, which indicates the number of rows of pixels. 
3. An integer, which indicates the number of columns of pixels. 
4. Eighty characters which hold an optional comment. 
The image data follows and is of type byte. [Ref: 4: page 15] 
2. Scanned Image file 
AS an image 1s scanned to identify on-road areas, this data along with the data 
used in identification of these areas is stored ina file. The same file format is used for 
both types of scanning data files. These files consist of an eighty-eight byte header of 
general data, followed by the data that identifies the on-road areas. Since the version 
of LISP used (see Chapter Four) allows only storage of integers one byte in length, the 
integers and floating point numbers are broken down into two bytes of data. Table 1 
shows the format of the first eighty-eight bytes of general data stored in the file. 
The next item stored is a single byte integer, the size of the window. The rest 
of the file contains the data on the identified on-road areas. Each identified on-road 
area has five bytes of data stored. The first byte is an integer indicating the selected 


feature or features utilized in determining that the area is on-road. The next two bytes 
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TABLE 1 
FILE FORMAT 


Field Name Bytes 


outer on road mean 

inner on road mean 

deviation of the outer on road mean 
deviation of the inner on road mean 
outer on road variance 

inner on road variance 

deviation of the outer on road variance 
deviation of the inner on road variance 
outer off road mean 

inner off road mean 

deviation of the outer off road mean 
deviation of the inner off road mean 
outer off road variance 

inner off road variance 

deviation of the outer off road variance 
deviation of the inner off road variance 
threshold of the outer mean 

threshold of the inner mean 

threshold of the outer variance 
threshold of the inner variance 
on_road variance product 

off road variance product 

on.road probability 

off road probability 

inner on road mean 

variance on the inner on road mean 
variance on the outer on road variance 
variance_on the inner on road variance 
inner off road mean 

variance on the inner off road mean 
variance of the outer off road variance 
variance of the inner off road variance 


hobo 


WWWWI WWW WI WW WWI GI WI HII WI DIDIWWDIPIW WDD WW 





are the x screen coordinate. The final two bytes are the y screen coordinate. These two 


byte integers are broken down as described above. 


C. SELECTED FEATURES 
1. Moments 

When using texture to identify image regions, there are two approaches to 
texture description: statistical and structural [Ref. 8: page 406]. The approach chosen 
for this work was statistical. In the statistical approach, there are many possible 
methods for extracting relative information about the textures. These methods include, 
but are not restricted to: means, weighted means, variance, weighted variances and 
more complete measurements such as inertia, local homogeneity and energy 


[Ref. 1,2: pages 368, 279]. The structural approach includes methods of pattern 
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matching using normalized correlation matching, and least-square error polynomial fits 
of image approximations [Ref. 1: page 369]. Since texture has meaning only when 
discussing a group of points or pixels, the image must be segmented into regions or 
areas. When discussing texture, these regions are referred to as windows. Again, there 
are several ways of dividing an image into windows. Windows can all be of uniform 
shape and size, or they can vary depending on prior analysis. In this study, all windows 
are square regions composed of between 4 and 225 pixels. Square windows were chosen 
for simplicity and the minimum and maximum size selected due to of the size of the 
roads on the maps used. The window that is being classified is defined as the inner 
window. The outer window refers to a region that is centered on the inner window, but 
is larger than the inner window. The outer window is used to help further classify the 
inner window. 

The three features chosen for this study were the mean of the inner window, 
the variance of the inner window and the variance of the outer window. The mean of 
the outer window was not used because early tests indicated the results would be 
virtually identical to the result from the inner window mean. The mean of the inner 
window defines the gray level of the window. The gray level is important because in 
most cases there is a distinct difference in this feature between the two textures. The 
variance from the inner window defines the granularity, or roughness of the window. 
The outer window is used to determine if the granularity varies as the size of the 
window approaches or exceeds the size of the object, in this case the road. These 
moments Were chosen because the calculations were relative simple and they represent 
features that intuitively suggest the different properties the human eye might pick out. 

2. Selecting Windows 

In order to detect on-road areas, the system must first be taught. The system 
is designed to distinguish between two type of textured areas. In order to instruct the 
system, the user must do several things. First, the user must decide on the size of the 
square inner and outer windows. The size of the windows depends on the size of the 
roads. The inner window should be smaller than the width of the roads and the outer 
window should be the of same size or larger than the width of the roads. The smaller 
window is the scan area being classified. When teaching the system, the inner window 
appears as a cursor to choose the on-road and off-road areas. 

The user must also determine the number of areas that are used to instruct the 
system. This number should be large enough to include all different types of on-road 


and off-road areas. 
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D. APPROACHES TO THRESHOLD SELECTION 
1. Selected Features Considered Separately 

As the image is scanned, one method to identify a window is to use each 
moment separately to determine the textures. In this first method, each scan window is 
tested against each of the three selected features, and results are stored indicating 
which of the features, if any, determine the scan window to be an on-road area. In 
order to test each feature, a threshold must be determined for each feature. Two 
methods to determine a threshold were tested. A simple method of picking the 
midpoint between on-road and off-road values was first used. This method produced 
good results. 

A second method to chose a threshold, which proved more successful and is 
the method now used in the system, uses the normal distribution and density functions 
to compute an optimal threshold. This equation. was derived as shown in Equations 
Per- 3.15. 


The Gaussian probability density function, also called the normal curve, is 





given by 
; —(x— 1)? 
2 
nc ee — No) 
N20. G (eqn 3.1) 


Where ft 1s the mean value of x and o is the standard deviation of x. This function is 
graphically depicted in Figures 3.1 and 3.2 

By definition, the shaded area under the left hand side of the curve in Figure 
3.1 is given by 


bk . (eqn 3.2) 
A= J Niu,o)dx 





Figure 3.1 Normal Density Function p(x). 


To determine A, it is necessary to standardize the Normal curve of N( Jt, o ) so that it 


is in the form N(0,1). To accomplish this, let 


y= 
Oo 
(eqn 3.3) 
Then, 
Z 
Z xX — 
yo — (eqn 3.4) 
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Substituting y~ into Equation 3.1 and this-result into Equation 3.2, it follows that 


ty ya 

1 2 

A= | —==e “ dy =erf(ty) 
! N20, (ecinesrs) 


where erf{x) is the error function or cumulative gaussian probability density function. 
Substituting the upper limit of Equation 3.2 into Equation 3.3, 





(eqn 3.6) 


In the same manner as above, the area A under the right side of the curve in Figure 3.1 
is given by 


A= [N(i1,0)dx = erf(t,) (eqn 3.7) 
ty’ ‘ 


where the lower limit t,- is defined by 


W=k+2(U-t,) =2u-t, (eqns; 3) 


Solving Equation 3.7 for t, 








t, = 2U — ty 
(eqn 3.9) 
Now substituting Equation 3.9 into Equation 3.6, 
Zz KH a H—ty 
Gg 
(eqn 3.10) 
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Now consider two distributions, N(My ©, ), and N(Hy 6» ), with a common 
threshold t, . Referring to Figure 3.2, the areas under the curves is given by 


A, =erf 





a 
3, (eqn 3.11) 





sae a 
02 


(eqn 3.12) 
If the probability of classifying an on-road pixel as off-road is the same as classifying 
an off-road pixel as on-road, then Ay = Eo) and thus, from Equation 3.11 and 3.12 


(Hy a5 ty) = (ty om [1>) 


Si, os (eqn 3.13) 


N( Hy To) 


N( Hy , 54 ) 





Figure 3.2. Normal Density Function of Two Gaussian Curves. 
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Equation 3.13 can be rewritten to obtain: 


(6, +0))t, = On, +0) 
1 2)tx Z 1 12 (eqn 3.14) 


Thus, solving for ty, from Equation 3.14, the final result for equal probability of these 


two types of error is 


_ OyPy + O2Hy 


0,+90 
oe (eames 5) 


[Ref. 9: pages 188-189}. In what follows, Equation 3.15 is the equation used to 
determine the individual thresholds of each of the three selected features. 
2. Selected Features Considered Together 
In this method of identifying a window, called a Gaussian scan, all three 
selected features are used together to determine if a region is on-road or off-road. In 
the methods described above, a classification based on a scalar measurement, a 
threshold was used for binary classification. In a more general classification based on 
vector measurements, a separating surface is needed instead of a simple threshold. 
There are two ways to approach this problem. The first uses a hyperplane and 
a measure of the relative importance of the two types of error. For example, if both 
types of error are equally important, then a plane is found that cuts equal volumes off 
both tails of the total distribution curves. A problem with this approach is that in 
general it doesn’t yield the minimum error. It yields the minimum error only in the 
special case of equal covariances for the Gaussian case [Ref. 10: page 194 Eq. 4]. 
_ A second approach that works in the more general case assumes that both 
types of error are equally important (1-0 loss function). In this approach,the optimal 
separating surface is a curved surface. The decision rule with respect to such a surface 


1s to select class i such that i maximizes the function, f{i), given by 


I 
. L 2 L (x.—m,.)? 
)= nF, Fh -= a ian 
. j=l oO 


jot ij (eGiieo mC) 


2] 


where, j = {inner window mean, outer window variance, inner window variance }, andi 
={ on-road, off-road }. The variable P; is the prior probability that a scan window is 
on-road or off-road. The second term in Equation 3.16 is called the variance product 
and is the product of the variances of the mean of the inner window, the variance of 
the outer window and the variance of the inner window. The last term in this equation 
is a generalized distance from the center of the class distributions. This function 
produces a minimal error rate for both types of misclassification errors. 

It should be noted that Equation 3.16 consists of all constants except for the 


X; Thus, the optimal choice has the shortest generalized distance from the class means. 


E. EXPERIMENTAL FACILITIES 
1. VAX/VMS System and the Digitizing Camera 
Connected to the VAX/VMS System is the EIKONIX 785 Digitizer Camera 
System. The camera is capable of digitizing an image up to 4096 X 4096 pixels and 
storing this image in a file. The camera is capable of taking either a color or black and 
white image. The camera is equipped with a 55mm SLR lens with Fstop range of 2.8 to 
22. The camera height is adjustable from 32 to 140 centimeters. 
2. Display and Scanning System 
The graphics system used is the Silicon Graphics Inc. IRIS Graphics 
Workstation. This system is capable of displaying images produced by the digitizing 
camera. The graphics station is limited to 1024 X 768 pixel display. Therefore, any 
image larger than this needs to be broken up. The images are displayed using RGB 


mode. 


F. SUMMARY 

This project was designed to find features that would aid in identifying on-road 
areas in aerial photographs. Two different methods of determining a threshold were 
examined. These methods were a simple binary method and a method where the 
probability of both types of possible errors were made equal. Also presented was a 
Gaussian method of combining several features to produce a minimum probability of 


Crror. 
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IV. IMPLEMENTATION OF THE RECOGNITION SYSTEM 


A. INTRODUCTION 

This chapter contains descriptions of the major LISP and C functions used in the 
road identification system. It includes descriptions of the function purpose, input 
variables, the results of the function call, and any restrictions inherent in each function. 


B. MAJOR SOFTWARE LANGUAGE AND MODULES 
1. Language 
Although Common LISP was the language preferred by the author for image 
processing, it was not available on a suitable host at the time this study was 
undertaken. Therefore, in what follows, the earlier dialect, Franz LISP (Ref. 11] will be 
used. Since Franz LISP is largely upward compatible with Common LISP, this is not 
viewed as a Serious limitation regarding further extensions of this study. 
2. LISP Functions 
a. function INIT-FILE 

The purpose of this function is to display the digitized map image on the 
IRIS graphics display. It requires as input the name of the file containing the digitized 
image. It also requires the x and y screen coordinates of the upper left corner where 
the image is to be displayed. Although this function displays an image up to 1024 X 
768 pixels, this system was designed to display two map images side by side. For this 
reason images should be no wider than 500 pixels. 

The result of calling this function is a displayed image of the digitized map. 
This function is normally called once with screen coordinates (0,766) when the next 
function to be called is main. This displays the map image in the correct location so 
that the instructions displayed on the graphics screen do not overlap the map. This 
function is called twice to display two identical maps side by side with screen 
coordinates (0,766) and (524,766) when the next function to be run is show-road. Show- 
road has an x coordinate input of either 0 or 524. 

The function init-file sets up the global variables lines and cols. Both of 
these values are read from the digitized map file. The C function display is the only 


function called from this function. 
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b. function MAIN 

The purpose of this function is to instruct the system in the difference 
between the two textures of on-road and off-road. The required inputs are the outer 
and inner window sizes and the number of points for the teaching sample. 

This LISP function is designed for interactive operator inputs. It displays a 
cursor the size of the inner window and requests that the operator choose the specified 
number of on-road and off-road areas. This function constructs a list of the sum of the 
points within each area, and from this list calculates the means and variances both on- 
road and off-road of the inner and outer windows. These values become global 
variables. This function also calculates the global variables of the variance product, 
used in the Gaussian function, the standard deviations and the thresholds for both the 
means and the variances. 

Other functions called in this function are the LISP functions 
a_window_mean, a_window_variance, mean, variance, and threshold. The C functions 
called are window!, and readsquare. 

¢. function SCAN-MAP 

The function scan-map scans the displayed image identifying on-road areas 
and stores these areas in two files. It is assumed that the map image is 500 X 500 pixels 
and is displayed at screen coordinates (0 766). The coordinates could be easily changed 
to input variables, but for the purposes of this study the use of the maximum square 
image was preferred. The inputs are the size of the outer and inner windows, the two 
file names for the output and the probability that any given area will be on-road. The 
two output files hold the results of the two methods used in identifying the scanned 
windows. The scan is conducted from left to right, top to bottom. If an area is 
identified as on-road, the lower right corner coordinate of the area is stored in the file, 
along with a code indicating which of the selected features was used to identify the 
area. 

The LISP functions called from this program are store-data, gaussian, and 
compare. The C functions are readsquare and writesquare. 

d. functions STORE-DATA and READ-DATA 
The functions store-data and read-data are complements of each other. The function 
store-data stores all the global variables in the file given as the input parameter. This 1s 
normally called from scan-map and is used to store the global data in the file before the 
scan data. The function read-data is used to read the data out of the file before 


displaying the scan results. It is normally called from show-road. 
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e. function SHOW-ROAD 
This function displays the results of the scan on the IRIS graphics 
workstation. The required inputs are the filename and the x screen coordinate at which 
the image is displayed. This function is designed to be called after init-file has been 
called to display the image. It assumes that the image was originally displayed and 


scanned with the y screen coordinate of 766. 


C. SUMMARY 

The functions described in this chapter were designed to be called in a specific 
order. However they can all be called independently, within the restrictions given and 
with the required global variables set. Once an image is scanned, and the results saved, 
the file can be displayed at any time, with or without superimposing it over the map 
image. If the system is taught on one map image, and there are several images to be 
scanned that the operator feels contain the same textures of the same size, the system 
need only be taught on the first map image, and this data can be used for all similar 


images, by first reading out of a file the global variables. 
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V. EVALUATION 


Aj. INTRODUCTION 

To evaluate the system, forty-eight different experiments were conducted. These 
experiments were broken down into four groups, depending on how the off-road 
regions were selected when teaching the system. Table 2 shows the breakdown of 
regions. Each group was broken down into two subgroups, a and b. Subgroup a was 
conducted with an outer window size of 10 X 10 and an inner window size of 5 X 5. 
Subgroup b had an outer window size of 9 X 9 and an inner window size of 3 X 3. 

All on-road regions chosen to instruct the system came from representative 
portions of both major roads shown in Figure 5.1 Specifically, the major roads are the 
vertical road and the road at approximately 20° from the horizonal. 


TAREE 
GROUPS FOR TEACHING OFF-ROAD REGIONS 


GROUP TYPE REGIONS 


15 verv d 


ark 
> very dark, 4 dark, 6 light, 2 very light 
ont 
15 very light 





B. 7 COLOR SCAN 

In this scan, all three selected features inner window mean, inner window 
variance, and outer window variance were individually tested. The threshold was 
chosen for each feature as shown in Equation 3.15 . The results were color coded (see 
Table 3) and displayed. 

The results from the first three groups of off-road regions, and subgroups from 
both sizes of windows were similar. The results can be seen in Figures 5.2, 5.3 and 5.4. 
These figures show that the major roads were found, with red indicating all three 
features were used. The last test, as seen in Figure 5.5, produced different results. When 


the larger windows were used, the inner variance picked out the road, the outer 
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TABLE 3 
7 COLOR SCAN COLOR CODES 


inner mean 

inner variance 

outer variance 

inner mean and inner variance 

inner variance and outer variance 

inner mean and outer variance 

inner mean outer variance and inner variance 





variance picked out the off-road and the mean picked out the entire image. Using the 
smaller windows, the mean picked out the road, while both variances picked out the 
off-road. In this map image, this result indicates that when choosing off-road regions 
to be very light, the difference in the mean for the larger windows was reduced to near 
zero. With the smaller windows, the mean was the only feature that gave good results. 
It is clear that the way in which off-road regions were taught, and to a lesser extent, 


the size of the windows, has a great deal of influence on the results in this type of scan. 


C. GAUSSIAN SCAN 

The Gaussian Scan is accomplished using all three selected features together as 
shown in Equation 3.16 . Forty experiments were run, using five different probabilities 
for a region being on-road. The four different groups given in Table 2 were used, as 
were the two different size of windows. Five gaussian scans were run in each subgroup 
using a priori road probabilities of .1, .3, .5,.7 and .9. » 

The results showed very little difference when the prior probabilities were 
changed within each subgroup. Each increase in probability added a few more windows 
identified as on-road. The difference in performance between the two sizes of windows 
was not significant except in the first test as seen in Figure 5.6. The smaller window 
sizes picked up a large amount of lighter off-road regions as being on-road. Figures 5.7, 
5.8 and 5.9 show varying degrees of results with the second and third test group 
showing the best results. 


D. CONCLUSION 

The procedure in running the experiments was to teach the system once for each 
subgroup and then to run the 7 color scan and the five Gaussian scans. As a result, 
the exact same off-road regions were used in each of these six scans. This made the 


data easier to compare. 


ai 


The results clearly show that the Gaussian type scan produced better results. In 
three of the four groups selected, the results were very good. The results also indicate 
that in the teaching portion, the regions chosen can improve the results. In this image, 
on-road areas are light, with some variance, and off-road regions vary from a few very 
light areas with no variance to very dark regions with no variance. It 1s important for 
the operator to choose either a varying sample of off-road regions as in group 2, or as 
in group 3, off-road areas that are close but not identical to on-road areas, for the best 


result. 
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VI. SUMMARY AND RECOMMENDATIONS FOR FURTHER 
RESEARCH 


A. SUMMARY 

This study was aimed at the practical problem of locating roads in aerial 
photographs. Rapid processing and accurate identification of road areas were the main 
objectives. To achieve these goals, simple texture features were chosen to identify the 
road areas. 

The first approach was to try a binary feature vector. In this approach, each of 
the three selected features was considered independently. The results of this approach 
were interesting, but not as good as hoped. The results indicated that the features were 
informative, but should be evaluated in a different way. 

A second approach was chosen. This approach was based on a more complex 
method of using all three selected features in combination in which the probability of 
misclassification is minimized. The results of this approach greatly improved the 
identification of road areas. 

Overall, the results obtained are useful in flagging areas for further attention by a 
human photo interpreter. However, window classification of pixels alone does not 
adequately separate road from non-road areas. Further processing using syntactic 


information is needed to obtain better classification of areas. 


B. RECOMMENDATIONS 

This work involved simple methods of classifying windows. More work on 
window classification is needed. A more complex method might be justified if the 
effectiveness of identifying road areas improved sufficiently. With anticipated 
improvements to the hardware of the Graphics Laboratory at NPS, it may be possible 
to use these more complex methods without a substantial increase in the processing 
time. 

The syntactic identification of road areas using an artificial intelligence expert 
system such as KEE or ART should be investigated. The system illustrated here could 
be used as a filter to locate possible road areas for the expert system to look at and 


further classify using syntactical features of road, such as shape and size. 
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A final recommendation is to investigate the possibility of improving human 
interaction in the identification of road areas. For example, a human could have 
available a low resolution aerial photograph of a large area. The computer would have 
the same area, but at a much higher resolution. This high resolution image could be 
processed by the computer. The computer could then flag areas for, the human photo 
expert to consider. The human could then confirm the computers identification using 
the low resolution photo, or could call up and display the area in question from the 
high resolution image. 
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APPENDIX A 
USER’S MANUAL 


l INTRODUCTION 

This document describes how to run the software to teach the system the 
difference between on-road and off-road areas, the software to scan the image to find 
the road, and the software to display the road areas. It also describes the file format of 
the digitized image. 


2, SOFTWARE DOCUMENTATION 
a. File Format 
The files containing the digitized map images are created by the EIKONIX 
camera and associated software on the VAX/VMS system and transferred to the IRIS 
(Ref. 4: page 15]. This application uses the files for black and white images only. Each 
file contains an eighty-six byte header followed by the image data. The header is 
constructed as follows: 


1. The first two bytes represent the file type. Format is integer. This will be 1’ for 
all black and white images. 


2. The third and forth bytes are the number of lines of image data. Format 1s 
integer. 


3. The fifth and sixth bytes are the number of columns of image data. Format 1s 
integer. 


4. The seventh through eighty-six byte are comments. Format 1s character. 
All integer fields are stored with the Least Significant byte(LSB).first. This application 
assumes all images are 500 columns wide. The image data immediately follows the 
header and consist of single bytes. Format is byte. The data is stored sequentially by 
line. For the file format of the scan files on the IRIS, see Table 1. 
b. Detailed User Instructions 

The software is written in Franz LISP and C. It is designed to be run from 
the IRIS side terminal. The display of the image and color coded on-road areas are 
directed to the associated IRIS display terminal. In order to run the program, all C 
functions must be compiled. It is recommended that all LISP functions be complied to 


Increase execution speed. 
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To run any of the software, it is necessary to be in the LISP environment, and 
to have the C functions and LISP function loaded. See Figure A.1 for an example of 


how to do this. 


(load ‘myfilename) 
(cfasl “display.o”. ”_display” ‘display_1 “c-function” ”-lgl”) 
(cfasl “ginit.o” ”_gint” ‘gint_| “c-function” ”-lgl”) 

(cfasl “windowl.o” ”_windowl” ’window1_1 ”c-function” ”-1gl”) 
(cfasl “readsquare.o” ” 


on 7 


_readsquare” ‘readsquare_| “c-function” ”-lgl”) 
(cfasl “writesquare.o” ”_writesquare” ‘writesquare_l ’c-function” ”-lgl”) 


oo 


(cfasl “colorcode.o” ”_colorcode” ‘colorcode_| “c-function” ”-lgl”) 





Figure A.l Sample of Loading Calls. 


I. Map Display 
To display the map image, call the LISP function init-file. The parameters 
are the filename of the file containing the digitized image, the x screen coordinate of 
the left side of the image, and the y screen coordinate of the top of the image. 
Example: (init-file ‘map.dig2 0 766). This will display map.dig2 in the upper right 
corner of the IRIS display terminal screen. 
2. Instruct Computer 
To teach the computer the difference between on-road and off-road areas 
on a given image, first display the image as in the example above, insuring that it is 
displayed at screen coordinates 0,766. Next call the LISP function main. Parameters 
are the inner window size, inner-window-size and number of points to teach with. 
Example: (main 9 3 15). It will then be necessary to use the left mouse button to 
select the points the computer request. Note the size of the selection box is the size of 
the outer-window:size. 
3. Scan Map 
To scan the map image, it is necessary to first generate the data produced 
when the system is taught. This can be accomplished in one of two ways. First of all 
the scan can be run immediately after the computer is taught, since the required data is 


assumed to be global. The second choice is to have saved the data from the teaching 
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phase in a file and retrieve that data. There are two LISP functions read-data and 
store-data that will accomplish this. Both functions have only the file name as a 
parameter. Once the data is global, the scanning phase is accomplished by running 
scan-map. Example: (scan-map outer-window-size inner-window-size outfilel outfile2 
prob) Outfilel is an output file which will contain the results of the seven color scan. 
Outfile2 is an output file which will contain the results of the gaussian scan. Prob is the 
probability the scanned window will be an on-road area. Note that the results of the 
gaussian scan will be displayed in the upper right of the screen. If desired, the image of 
the map can be displayed there before running the program. This will result in an 
overlay effect. To display the map use the following: (init-file filename 510 766) 
4, Display Road 

To display the road the LISP function show-road is used. Parameters are 
the filename of a file created in scam-map and the offset. The road can only be 
displayed in the upper part of the screen. The offset is the x screen coordinate. 


Example: (show-road map.gaussian 510). 
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APPENDIX B 
PROGRAM 


This set of functions is designed to identify on-road 


areas in digitized arial photographs 


GLOBAL VARIABLES: 
outer_on_road_mean 
inner_on_road_mean 
outer_on_road_var 
inner_on_road_var 
outer_off_road_mean 
inner_off_raod_mean 
outer_off_road_var 


inner_off_road_var 


outer_var_on_mean 
inner_var_on_mean 
outer_var_on_var 
inner_var_on_var 
outer_var_off_mean 
inner_var_off_mean 
outer_var_off_var 


inner_var_off_var 


var_prod_on_road : 


var_prod_off_road : 


mean of 


mean of 


mean of 


mean of 


mean of 


mean of 


mean of 


mean of 


variance 
variance 
variance 
variance 
variance 
variance 
variance 


variance 


on-road means of outer windows 
on-road means of inner windows 
on-road variances of outer windows 
on-road variances of inner windows 
off-road means of outer windows 
off-road means of inner windows 
off-road variances of outer windows 


off-road variances of inner windows 


of on-road means of outer windows 
of on-road means of inner windows 
of on-road variances of outer windows 
of on-road variances of inner windows 
of off-road means of outer windows 
of off-road means of inner windows 
of off-road variances of outer windows 


of off-road variances of inner windows 


variance product is the ln of the square root 


of inner_var_on_mean * outer_var_on_var * 


inner_var_on_var 


variance product is the ln of the square root 


of inner_var_off_mean * outer_var_off_var * 


inner_var_off_var 
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outer_dev_on_mean : 


inner_dev_on_mean 


outer_dev_on_var =: 
inner_dev_on_var :; 
outer_dev_off_mean: 
inner_dev_off_mean: 
outer_dev_off_var : 


inner_dev_off_var : 


outer_mean_threshold 
inner_mean_threshold 
outer_var_threshold 


inner_var_threshold 


standard 
standard 
standard 
standard 
standard 
standard 
standard 


standard 


dev 
dev 
dev 
dev 
dev 
dev 
dev 


dev 


: threshold 
: threshold 


: threshold 


: threshold 


of 
of 
Of 
of 
of 
of 
of 
of 


mean of on-road outer windows 

mean of on-road inner windows 
variance of on-road outer windows 
variance of on-road inner windows 
mean of off-road outer windows 
mean of off-road inner windows 
variance of off-road outer windows 


variance of off-road inner windows 


for mean of outer windows 


for mean of inner windows 


for variance of outer windows 


for variance of inner windows 
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o AAKKAAAKKAAKAKKKAKKKAAAK KKK AKKAKAKKRKAKKRKAAKKAAKKAAAKK AK AKKK ARK 


; MAIN * 
; Function is to instruct system in the differance between * 
; on-road and off-road areas * 
; Inputs: * 
* outer_window_size: size of outer scan window * 
; inner_wondow_size: size of inner scan window * 
; points : number of points used to teach * 
; Functions: * 
. window_l * 
; readsquare_l * 
; a_window_mean zx 
; a_window_variance * 
; mean * 
; variance x 
; threshold * 


e RRKAARKKAAAKKKKKKKKAAKKK KKK AKRKKKAKRKKKKAKKRKKKAKAAAKKKAAKKRKKKAARK 


(defun main (outer_window_size inner_window_size points) 
(setq outer_cursorval (new-vectori-word (expt outer_window_size 


2))) 


(setq inner_cursorval (new-vectori-word (expt inner_window_size 


2))) 


(setq xcursor (new-vectori-word 2)) 


(setq ycursor (new-vectori-word 2)) 


ax 

;***X Loop to learn on road areas 

ok 

mcon<( 1 0 (+ i 1))) S 

((= i points) t) 

i 

;*xAAKK read in pixel values for the on road windows 

i 
(windowl_l outer_window_size lines cols xcursor ycursor 1) 


(readsquare_l (vrefi-word xcursor 0) (vrefi-word ycursor 0) 


45 


outer_window_size outer_cursorval) 

(vseti-word xcursor 1(+ (+(-(/ outer_window_size 2) 

(/ inner_window_size 2))1) 
(vrefi-word xcursor 0))) 
(vseti-word ycursor 1(+ (+(-(/ outer_window_size 2) 

(/ inner_window_size 2))1) 
(vrefi-word ycursor 0))) 
(readsquare_l (vrefi-word xcursor 1) (vrefi-word ycursor 1) 

inner_window_size inner_cursorval) 
(cond ((= i 0) 

(setq outer_on_mean_list (list (a_window_mean 

outer_window_size outer_cursorval))) 

(setq inner_on_mean_list (list (a_window_mean 

inner_window_size inner_cursorval) ) ) 

(setq outer_on_var_list (list (a_window_variance 
outer_window_size 
outer_cursorval 
(car outer_on_mean_list) 

oy 

(setq inner_on_var_list (list (a_window_variance 
inner_window_size 
inner_cursorval 
(car inner_on_mean_list) 


)) =) 
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((> i 0) 
(attach (a_window_mean outer_window_size 
outer_cursorval) outer_on_mean_list) 
(attach (a_window_mean inner_window_size 
inner_cursorval)inner_on_mean_list) 
(attach (a_window_variance outer_window_size 
outer_cursorval (car outer_on_mean_list) ) 
outer_on_var_list) 
(attach (a_window_variance inner_window_size 
inner_cursorval (car inner_on_mean_list) ) 


inner_on_var_list)  ) 


ox 


a 


;xkXKX Loop to learn off road areas 
ox 


, 


feo (( 1 0 (+ i 1))) 
((= i points) t) 
a 
;xxAKX read in pixel values for the off road windows 


ox 


? 


(windowl_1l outer_window_size lines cols xcursor ycursor 2) 
(readsquare_l (vrefi-word xcursor 0) (vrefi-word ycursor 0) 
outer_window_size outer_cursorval) 

(vseti-word xcursor 1(+ (+(-(/ outer_window_size 2) 

(/ inner_window_size 2))1) 
(vrefi-word xcursor 0))) 
(vseti-word ycursor 1(+ (+(-(/ outer_window_size 2) 

(/ inner_window_size 2))1) 
(vrefi-word ycursor 0))) 
(readsquare_l (vrefi-word xcursor 1) (vrefi-word ycursor 1) 


inner_window_size inner_cursorval) 
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(cond ((= i 0) 

(setq outer_off_mean_list (list (a_window_mean 

outer_window_size outer_cursorval))) 

(setq inner_off_mean_list (list (a_window_mean 

inner_window_size inner_cursorval) )) 

(setq outer_off_var_list (list (a_window_variance 
outer_window_size 
outer_cursorval 

(car 
outer_off_mean_list) 
yy 

(setq inner_off_var_list (list (a_window_variance 

inner_window_size 

inner_cursorval 

(car 
inner_off_mean_list) 


))) ) 
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(rio) 

(attach (a_window_mean outer_window_size 
outer_cursorval) outer_off_mean_list) 

(attach (a_window_mean inner_window_size 
inner_cursorval) inner_off_mean_list) 

(attach (a_window_variance outer_window_size 
outer_cursorval (car outer_off_mean_list) ) 
outer_off_var_list) 

(attach (a_window_variance inner_window_size 
inner_cursorval (car inner_off_mean_list) ) 


inner_off_var_list)  ) 


ox 

,xxxxKK Figure mean of given list 

+ 

(setq outer_on_road_mean (mean outer_on_mean_list points) ) 
(setq inner_on_road_mean (mean inner_on_mean_list points) ) 
(setq outer_on_road_var (mean outer_on_var_list points) ) 
(setq inner_on_road_var (mean inner_on_var_list points) ) 
(setq outer_off_road_mean (mean outer_off_mean_list points) ) 
(setq inner_off_road_mean (mean inner_off-mean_list points) ) 
(setq outer_off_road_var (mean outer_off_var_list points)) 
(setq inner_off_road_var (mean inner_off_var_list points) ) 


ox 


7 


;xxAKK Figure variance of given list 
ox 


4 


(setq outer_var_on_mean (variance outer_on_mean_list points 
outer_on_road_mean) ) 
(setq inner_var_on_mean (variance inner_on_mean_list points 


inner_on_road_mean) ) 


(setq outer_var_on_var (variance outer_on_var_list points 


outer_on_road_var)) 
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(setq inner_var_on_var (variance inner_on_var_list points 
inner_on_road_var) ) 
(setq outer_var_off_mean (variance outer_off_mean_list points 
outer_off_road_mean) ) 
(setq inner_var_off_mean (variance inner_off_mean_list points 
inner_off_road_mean) ) 
(setq outer_var_off_var (variance outer_off_var_list points 
outer_off_road_var) ) 
(setq inner_var_off_var (variance inner_off_var_list points 


inner_off_road_var) ) 


oX* 


c 


;xxkKK Figure the variance product to be used 


;xkAKK in the gaussian function : 
ok 
e 


(setq var_prod_on_road (log (sqrt (float (product inner_var_on_mean 


outer_var_on_var inner_var_on_var))))) 


(setq var_prod_off_road (log (sqrt (float (product inner_var_off_mean 


outer_var_off_var inner_var_off_var))))) 
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xX 


s 


,kAAKK figure standard dev of given list 


7% 


, 


(setq 
(setq 
(setq 
(setq 


(setq 
(setq 
(setq 
(setq 


ox 


c 


outer_dev_on_mean (sqrt (float 
inner_dev_on_mean (sqrt (float 
(sqrt (float 
(sqrt (float 


outer_dev_on_var 


inner_dev_on_var 


outer_dev_off_var 


inner_dev_off_var 


outer_var_on_mean) ) ) 
inner_var_on_mean) ) ) 
outer_var_on_var) )) 


inner_var_on_var) )) 


outer_dev_off_mean (sqrt (float outer_var_off_mean) ) ) 
inner_dev_off_mean (sqrt (float inner_var_off_mean) ) ) 
(sqrt (float outer_var_off_var) )) 


(sqrt (float inner_var_off_var) )) 


,xxkKK Figure thresholds for means and vars 


7% 


c 


(setq outer_mean_threshold (threshold outer_dev_off_mean 


outer_dev_on_mean 
outer_on_road_mean 


outer_off_road_mean) ) 


(setq inner_mean_threshold (threshold inner_dev_off_mean 


inner_dev_on_mean 
inner_on_road_mean 


inner_off_road_mean) ) 


(setq outer_var_threshold (threshold outer_dev_off_var 


outer_dev_on_var 
outer_on_road_var 


outer_off_road_var) ) 


(setq inner_var_threshold (threshold inner_dev_off_var 


inner_dev_on_var 
inner_on_road_var 


inner_off_road_var) ) 


g RAEKKAKAKKAKKKAAKAKKKKKAKKKKKKAAKKKKKAKKKRKKKAKAKKKKAAKAKKKKAKAKRKKAKAKK 


; INIT-FILE 

; Displays image 

; Inputs: 

; filename: name of file holding digitized image 

; Gisplay-at-x: screen x coordinate (left side of image) 
; display-at-y: screen y coordinate (top) 


; Functions: 


e He He Fe ee eH 


; display_l 
5 KAKKAAAKKAKAARAKRARAKKAAKAKEKKKARAK RARKKRAK KARA KKK RRA RKR KK 


(defun init-file (filename display-at-x display-at-y) 
(display_l (get_pname filename) display-at-x display-at-y) 
(setq fd (infile filename) ) 
(fseek fd 2 0) 
(setq lines (+ (tyi fd) (* 256) (tvierad 
(setq cols (+ (tyi fd) (* 256°(tv192a) 
(close fd) 
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; A_WINDOW_MEAN * 
; Finds the mean of a given window x 
; Inputs: * 
> Windowsize: the size of the window x 
; cursorval: the array holding the values of each pixel * 
; in the window k 
; Returns: * 

z 


; the mean of the given window 
5p RAKAKAAAKAAKAARKAAKKARKARKKARKAARKAKKKK KKK KAKA KKK RK RK KK 


(defun a_window_mean (windowsize cursorvalL) 
mado ({ j 1 (+ j 1)) 
(sum (vrefi-word cursorval 0) 
(+ sum (vrefi-word cursorval j)))) 
((= j (expt windowsize 2)) sum) 
)(expt windowsize 2)) 
) 


s RAKRARKRKAKAKRKARKARKARARKRAKRKRKARKARAKARRAK AK AKA RRKK RRA AKARRKAKEK 


; A_WINDOW_VARIANCE * 
; find variance in a window * 
; Inputs: * 
; Wwindowsize: the size of the given window * 
; cursorval: an array holding the value of each pixel * 
; in the window * 
; Window_mean: the mean of the given window x 
; Returns: * 
; the variance of the given window * 


g KAAKKAAKKAKKKKKKKKKKKKKKAAKAKKKKKAKKKKKAKKKAKKKAKKKKAKKKAKKKKKKKKKK 


(defun a_window_variance (windowsize cursorval window_mean) 
mm ao ((7 1 (+ 5 1)) 
(sumsq (expt(vrefi-word cursorval 0) 2) 
(+ sumsq (expt(vrefi-word cursorval j) 2)))) 
((= j (expt windowsize 2)) sumsq) 


)(expt windowsize 2))(expt window_mean 2)) 
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, KAKKAARAKAKKRRKKAAAKR KKK RAK RAK KKARAKRKKAKAARAKRAKARRARRAR ARK EK 
: MEAN 

; find mean of the input list 
; Inputs: 


list: a list containing the values you wish to find the 


=e 


; the mean of 
; points: the number of points in the list 


» Returns: 


i, a, a, a, i, a, a, » 


; the mean of the given list of points 
s KAKKKKKKRARKRKRKRK KARA RARER RKRKIKAK RAR ARAR ARK RAR AAA 
(defun mean (list points) 
(/ ( do (( sum (car list) (+ sum (car rest_of_list))) 
(rest_of_list (cdr list) (cdr restsofmmce 
(( null rest_of_list) sum) ) points) 


) 
, RAERKAAKAKAKKAKKKRRKRRAKRKKRRARRAKA RAKE RAK RAK RRA RK ARRKERER 
> VARIANCE x 
; find variance of given list and the mean of * 
; the list, this is not normalized * 
; Inputs: x 
; aS : a list of points you wish to find the x 
: variance of = 
; points : the number on points in list * 
; the_mean: the mean of the list x 
> Returns; * 
“ 


; the variance of a list of points 

p KREAKRKAAKARARKAK RAK AK KAKA AKRAK RAR RAK KK RAR RRR KARA RAR REI 

(defun variance (list points the_mean) 
(quotient (do ((sum (expt (diff (car list) the_mean) 2) 

(add sum (expt (diff(car rest_of_list) 

the_mean) 2))) 
(rest_of_list (cdr list) (cdr rest of ivsee 

((null rest_of_list) sum) ) points ) 
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g RARKKKKK KK KA KRKKAKAKK RKRKKKAKAKKKKKAKKKKAKAK AKA RK KAAKKKAKKKK RK 


; Returns: 


; THRESHOLD * 
; find the threshold given two means and two standard dev x 
; Inputs: * 
: dev_off: standard dev of off road area 3 
; dev_on : standard dev of on road area x 
; OLE : mean of off road area x 
; on s mean of on road area x 

* 

* 


; the threshold between the on-road and off-road 
4g RARKRAKKKAKRAAKK AK KAKA AKER RRK KAR AK RR KAR AKA RK KRAKKR RAK KR KAK 


(defun threshold (dev_off dev_on off on) 
(quotient (add (product dev_off on) 
(product dev_on off)) 
(add dev_on dev_off) ) 


0 


g RARKKKKAKKAKAAKKKKRKKAKAAKRKKAKAKAKAAKKKKAKKKKKK AK AKKKKKRRKKKARKRKKKARKRE 


; COMPARE x 
> Given inner or outer, on and off road means or variances * 
; and the threshold returnes true if on road area x 
; Inputs: x 
; Window ; inner or outer, mean or var, for a window x 


; threshold : threshold for the given window* 
; true_type : inner or outer, mean or var, for on road* 


; false_type: inner or outer, mean or var, for off road x 
5 KAKKARKARKAK AK KRKK KIKI R RARER RK AK AK AKA RRRIK IRE 


(defun compare (window threshold true_type false_type) 
(cond (( or (and (> window threshold) 
(> true_type false_type) 
) 
(and (< window threshold) 
(< true_type false_type) 


(t 
() 
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 RAKKKAKKAKKKKAKKKKKKKKKKAKKKKAKKKKKKKKKKKAKKKKAKKAKAKAAKAKK KK 


; GAUSSIAN ~ 


finds the gaussian property of either on-road or off-road * 


™=e 


; areas. 
; Inputs: * 
; prob : probality an area is on-road * 
; var-prod . : variance product * 
; window_inner_mean: mean of the given inner window x 
: window_outer_var : var of the given outer window x 
; window_inner_var : var of the giver inner window x 
; mean_inner_mean : mean of all the inner windows means * 
; mean_outer_var s mean of all the outer window vars * 
. mean_inner_var s mean of all the inner window vars * 
; var_inner_mean : var of all the inner window means * 
; var_outer_var : var of all the outer window vars * 
; var_inner_var ; var of all the inner window vars! * 
> Returns: * 
; the gaussian value * 


g RAKARKKARKKKKKKKKKKKAKKKKKKKK KKK AKKKKAKKKKKRKKAKRKKKRKKKAKKKKK KKK 


(defun gaussian (prob var_prod 
window_inner_mean window_outer_var 
window_inner_var mean_inner_mean 
mean_outer_var mean_inner_var 
var_inner_mean  var_outer_var  var_inner_var) 
(diff (diff (log prob) var_prod) 
(quotient (add (quotient (expt (diff window_inner_mean 
mean_inner_mean) 2) var_inner_mean) 
(quotient (expt (diff window_outer_var 
mean_outer_var ) 2) var_outer_var ) 
(quotient (expt (diff window_inner_var 


mean_inner_var ) 2) var_inner_var ) )2)) 


a7 


s RRKAKKAAKKARKKAKKKAKKA KKK AKER KA KRKRARKRAKKKKKARKAKKKKRKKARAAKKRKR KK 


; SCAN-MAP * 
; scans the displayed image, identifies on-road areas and * 
; stores these areas ina file. * 
; Inputs: * 
; outer_size : the size of the outer window x 
; inner_size ‘3; the size of the inner window * 
; outfilenamel : file to hold on-road areas from 7 color * 
; scan * 
> outfilename2 : file to hold on-road areas from gaussian * 
; scan x 
7 Dros : probality an area is on-road * 
; Functions: * 
; readsquare_l * 
; gaussian * 
; compare * 
; Wwritesquare_l * 
> store-data * 


g RARKRKAAAKKAKKAAAKRKKRAKAKRKKKAA RK RKAAAA RK RAK AK RK RAK RRA RRA KAA 


(defun scan-map (outer_size inner_size outfilenamel outfilename2 
prob) 
(setq fdl (outfile outfilename1) ) 
(setq fd2 (outfile outfilename2) ) 
(setq on_road_prob prob) 
(store-data fd1) 
(store-data fd2) 
(tyo inner_size fd1l) 
(tyo inner_size fd2) 
(setq outer_windowvector(new-vectori-word (expt outer_size 2))) 
(setq inner_windowvector(new-vectori-word (expt inner_size 2))) 
(do (( i (+ (- 766 outer_size ) 1) (- i inner_size))) 
((< i (+ (- 766 lines) 1)) t) 
(setq inner_i (+ i (+ (- (/ outer_size 2)(/ inner_size 2)) 
1))) 


(msg "scan is now scanning map line " i N) 
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(do (( 3 0 (+ j inner_size))) 
((> j (- cols outer_size))t) 


(setq inner_j (+ j (+ (- (/ outer_size 2)(/ inner_size 2)) 
Do) 


(readsquare_l j i outer_size outer_windowvector) 
(readsquare_1l inner_j inner_i inner_size 


inner_windowvector ) 


(setq outer_window_mean (a_window_mean outer_size 
outer_windowvector) ) 
(setq inner_window_mean (a_window_mean inner_size 
inner_windowvector) ) 
(setq outer_window_var (a_window_variance outer_size 
outer_windowvector 


outer_window_mean) ) 


(setq inner_window_var (a_window_variance inner_size 
inner_windowvector 


inner_window_mean) ) 


Do 


0X 


, 


;kxkKKK the gaussian method for the scan 
ok 


f 


(cond ((> (gaussian on_road_prob var_prod_on_road 
inner_window_mean outer_window_var 
inner_window_var inner_on_road_mean 
outer_on_road_var inner_on_road_var 
inner_var_on_mean outer_var_on_var 
inner_var_on_var) 

(gaussian (diff 1 on_road_prob) var_prod_off_road 
inner_window_mean outer_window_var 
inner_window_var inner_off_road_mean 
outer_off_road_var inner_off_road_var 
inner_var_off_mean outer_var_off_var 
inner_var_off_var) 

) 

(tyo 7 £d2) 

(tyo (car (Divide inner_j 256)) fd2) 

(tyo inner_j fd2) 

(tyo (car (Divide inner_i 256)) fd2) 

(tyo inner_i fd2) 

(writesquare_l (+ 510 inner_j) inner_i inner_size 255 0 0) 


)) 
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* 


° 
a 


;**AXX scanning with seven color output 


i 


a 


(cond ( (compare inner_window_mean inner_mean_threshold 


inner_on_road_mean inner_off_road_mean) 


(cond ( (compare inner_window_var 


inner_on_road_var 


inner_var_threshold 


inner_off_road_var) 


(cond ( (compare outer_window_var 


outer_var_threshold 


outer_on_road_var 


outer_off_road_var) 


-(writesquare_l (+ 510 inner_j) inner_i inner_size 255 0 0) 


) 


GEvO 7 LdL) 

(tyo (car (Divide inner_j 256)) 
fd1) 

(tyo inner_j fdl) 

(tyo (car (Divide inner_i 256)) 
Pal) 

(tyo inner_i fdl) 


Cee 


;(writesquare_l (+ 510 inner_j) inner_i inner_size 0 255 255) 


) 


(tyo 4 fdl) 

(tyo (car (Divide inner_j 256)) 
fd1) 

(tyo inner_j fdl) 

(tyo (car (Divide inner_i 256)) 
lt) 

(tyo inner_i fd1) 


( (cond ( (compare outer_window_var 
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outer_var_threshold 
outer_on_road_var 
outer_off_road_var) 
;(writesquare_l (+ 510 inner_j) inner_i inner_size 150 100 100) 
(tyo Ge 1d) 
(tyo (car (Divide inner_j 256)) 
fda 
(tyo inner_j fdl) 
(tyo (car (Divide inner_i 256) ) 
fat) 
(tyo inner_i fd1l) 


) 
(Fe 
;(writesquare_l (+ 510 inner_j) inner_i inner_size 0 255 0) 
Cryo 1 fdl) 
(tyo (car (Divide inner_j 256)) 
fdl) 
(tyo inner_j fdl) 
(tyo (car (Divide inner_i 256)) 
fd1) 
(tyo inner_i fdl) 
) 


) 

( (cond ( (compare inner_window_var 
inner_var_threshold 
inner_on_road_var 
inner_off_road_var) 

(cond ( (compare outer_window_var 
outer_var_threshold 
outer_on_road_var 
outer_off_road_var) 


;(writesquare_l (+ 510 inner_j) inner_i inner_size 255 0 255) 
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) 


) 


ctyo 5 fd) 

Civomtcar (Davide inner} 256) ) 
fal ) 

(tyo inner_j fdl) 

(tyo (car (Divide inner_i 256)) 
pal) 

(tyo inner_i fdl) 


(Se 


;(writesquare_l (+ 510 inner_j) inner_i inner_size 0 0 255) 


(ayo. 2 £dly) 
(tyo (car (Divide inner_j 256)) 
Cay) 
(tyo inner_j fdl) 
(tyo (car (Divide inner_i 256)) 
fdl) 
(tyo inner_i fdl) 


( (compare outer_window_var 


outer_var_threshold 
outer_on_road_var 


outer_off_road_var) 


;(writesquare_l (+ 510 inner_j) inner_i inner_size 255 255 0) 


(tyo 
(tyo 
(tyo 
(tyo 
(tyo 


3 fal) 
(car (Divide inner_j 256)) fdl) 


inner_j fdl) 
(car (Divide inner_i 256)) fdl) 


inner_i fdl) 


) Dem etclose fdl) (close £d2) ) 
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s KARKKAAKAAAARAKAKKAKKAKKEKKRKKAAKRKARKARK RRA RKRARKRKKKERERKKKK 
; STORE-DATA 

; fixes the data up so it can be placed in a file. 
; This version of LISP only stores fixnumbers 


PLA UES - 


~ Ft Ft Fe + 


; fd: filename of file to store data in 
5 KKAKARAKAKARARKKRAKAKAKAK KKK KRKK KKK KR RK KAKA RRR RAK ERR RR IK 


(defun store-data (fd) 
(tyo (car (Divide (fix outer_on_road_mean) 256)) fd) 
(tyo (fix outer_on_road_mean) fd) 
(tyo (car (Divide (fix inner_on_road_mean) 256)) fd) 


(tyo (fix inner_on_road_mean) fd) 


(tyo (fix (quotient (product outer_dev_on_mean 1000) 65536)) fd ) 

(tyo (fix (quotient 
[diff (product outer_dev_on_mean 1000) (product 
(fix (quotient (product outer_dev_on_mean 1000) 
65536) ) 65536)-—"256)) "fa, 

(tyo (fix (product outer_dev_on_mean 1000)) fd) 


(tyo (fix (quotient (product inner_dev_on_mean 1000) 65536)) 
fd) 
(tyo (fix (quotient 
[diff (product inner_dev_on_mean 1000) (product 
(fix (quotient (product inner_dev_on_mean 1000) 
65536)) 65536)- 256)) iar 
(tyo (fix (product inner_dev_on_mean 1000)) fd) 


(tyo (car (Divide (fix outer_on_road_var) 256)) fd) 


(tyo (fix outer_on_road_var) fd) 


(tyo (car (Divide (fix inner_on_road_var) 256)) fd) 


(tyo (fix inner_on_road_var) fd) 


(tyo (fix (quotient (product outer_dev_on_var 1000) 65536)) fd) 
(tyo (fix (quotient 
[diff (product outer_dev_on_var 1000) (product 
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(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 
(tyo 
(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 
(tyo 


(tyo 
(tyo 


(fix (quotient (product outer_dev_on_var 1000) 
65536)) 65536)- 256)) fd) 
(fix (product outer_dev_on_var 1000)) fd) 


(fix (quotient (product inner_dev_on_var 1000) 65536)) fd) 
(fix (quotient 
[diff (product inner_dev_on_var 1000) (product 
(fix (quotient (product inner_dev_on_var 1000) 
65536)) 65536)- 256)) fd) 
(fix (product inner_dev_on_var 1000)) fd) 


(car (Divide (fix outer_off_road_mean) 256)) fd) 
(fix outer_off_road_mean) fd) 
(car (Divide (fix inner_off_road_mean) 256)) fd) 


(fix inner_off_road_mean) fd) 


(fix (quotient (product outer_dev_off_mean 1000) 65536)) fd) 
(fix (quotient 
[diff (product outer_dev_off_mean 1000) (product 
(fix (quotient (product outer_dev_off_mean 1000) 
65536)) 65536)- 256)) fd) 
(fix (product outer_dev_off_mean 1000)) fd) 


(fix (quotient (product inner_dev_off_mean 1000) 65536)) fd) 
(fix (quotient 
[diff (product inner_dev_off_mean 1000) (product 
(fix (quotient (product inner_dev_off_mean 1000) 
G5 SGC py moo saG)— 9 256)) itd) 
(fix (product inner_dev_off_mean 1000)) fd) 


(car (Divide (fix outer_off_road_var) 256)) fd) 


(fix outer_off_road_var) fd) 


(car (Divide (fix inner_off_road_var) 256)) fd) 


(fix inner_off_road_var) fd) 


(fix (quotient (product outer_dev_off_var 1000) 65536)) fd) 
(fix (quotient 
[diff (product outer_dev_off_var 1000) (product 
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(tyo 


(tyo 
(tyo 


(tyo 


(tyo 


(tyo 


(tyo 


(tyo 


(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(fix 


(fix (quotient (product outer_dev_off_var 1000) 
65536)) 65536)- 256)) £d) 
(product outer_dev_off_var 1000)) fd) 


(fix (quotient (product inner_dev_off_var 1000) 65536)) fd) 


Cra 


CEix 


(fix 


(fix 


(fix 


(fix 


(fix 


(fix 


(quotient 
_ (diff (product inner_dev_off_var 1000) (product 
(fix (quotient (product inner_dev_off_var 1000) 
65536)) 65536)- 256)) fd) 
(product inner_dev_off_var 1000)) fd) 


(quotient (product outer_mean_threshold 1000) 
65536)) fd) 
(quotient 
[diff (product outer_mean_threshold 1000) (product 
(fix (quotient (product outer_mean_threshold 
1000)65536)) 65536)- 256)) fd) 
(product outer_mean_threshold 1000)) fd) 


(quotient (product inner_mean_threshold 1000) 
65536), mec) 
(quotient 
[diff (product inner_mean_threshold 1000) (product 
(fix (quotient (product inner_mean_threshold 
1000)65536)) 65536)- 256)) fd) 
(product inner_mean_threshold 1000)) fd) 


(fix (quotient (product outer_var_threshold 1000) 65536)) fd) 


(fix 


(tas 


(quotient 
[diff (product outer_var_threshold 1000) (product 
(fix (quotient (product outer_var_threshold 
1000)65536)) 65536)] 256)) fd) 
(product outer_var_threshold 1000)) fd) 


(fix (quotient (product inner_var_threshold 1000) 65536)) fd) 


(fix 


(quotient 
[diff (product inner_var_threshold 1000) (product 


(fix (quotient (product inner_var_threshold 
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(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 
(tyo 


(tyo 


(tyo 


(tyo 
(tyo 


1000)65536)) 65536)=- 256)) fd) 


(fix (product inner_var_threshold 1000)) fd) 


(fix (quotient (product var_prod_on_road 1000) 65536)) fd) 


(£1x 


(fix (product var_prod_on_road 1000)) fd) 


(quotient 


[diff (product var_prod_on_road 1000) (product 
(fix (quotient (product var_prod_on_road 1000) 


o5586)) 65536)= 256)) £d) 


(fix (quotient (product var_prod_off_road 1000) 65536)) fd) 


(fix 


(fix 


(£1x 
(fix 


(£1x 
(fix 


(£1x 


(quotient 


[diff (product var_prod_off_road 1000) (product 
(fix (quotient (product var_prod_off_road 1000) 
e5536) oo 5s6) | 256)) £d) 
(product var_prod_off_road 1000)) fd) 


(quotient (product on_road_prob 1000) 65536)) fd) 


(quotient 


[diff (product on_road_prob 1000) (product 
(fix (quotient (product on_road_prob 1000) 65536)) 


(product on_road_prob 1000)) fd) 


65536) ] 


(quotient (product (diff 1 on_road_prob) 1000) 


(quotient 


65536)) fd) 


[diff (product (diff 1 on_road_prob) 1000) 
(product (fix (quotient (product (diff 1 
on_road_prob)1000)65536)) 65536)] 256)) fd) 
(fix (product (diff 1 on_road_prob) 1000)) fd) 


256)) fd) 


(fix (quotient (product inner_on_road_mean 1000) 65536)) fd) 


(fix (quotient 


[diff (product inner_on_road_mean 1000) (product 


(fix (quotient (product inner_on_road_mean 1000) 


65536)) 65536) ] 


6/7 


256)) £d) 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 
(tyo 


(tyo 


(tyo 


(fix 


CPi x 
(fix 


(fix 


(fix 
(fix 


(fix 


(product inner_on_road_mean 1000)) fd) 


(quotient inner_var_on_mean 65536)) fd) 


(quotient 
(diff inner_var_on_mean (product 
(fix (quotient inner_var_on_mean 65536)) 65536)) 
51200) Ed) 
inner_var_on_mean) fd) 


(quotient outer_var_on_var 65536)) fd) 


(quotient 
(diff outer_var_on_var (product 
(fix (quotient outer_var_on_var 65536)) 65536)) 
256) £0) 
outer_var_on_var ) fd) 


(fix (quotient (product inner_var_on_var 1000) 65536)) fd) 


(£1x 


(fix 


(quotient 
[diff (product inner_var_on_var 1000) (product 
(fix (quotient (product inner_var_on_var 1000) 
65536)) 65536) }ezso)) wea 


(product inner_var_on_var 1000)) fd) 


(fix (quotient (product inner_off_road_mean 1000) 65536)) fd) 


(fix 


(f1x 


(fix 
(fix 


(fix 


(fix 


(quotient 
(diff (product inner_off_road_mean 1000) (product 
(fix (quotient (product inner_off_road_mean 1000) 
65536)) 65536)) 256))eeear 
(product inner_off_road_mean 1000)) fd) 


(quotient inner_var_off_mean 65536)) fd) 


(quotient 
(diff inner_var_off_mean (product 
(fix (quotient inner_var_off_mean 65536)) 65536)) 
256)) fd) 
inner_var_off_mean) fd) 


(quotient (product outer_var_off_var 1000) 65536)) 
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fd) 
(tyo (fix (quotient 
(diff (product outer_var_off_var 1000) (product 
(fix (quotient (product outer_var_off_var 1000) 
65536)) 65536)) 256)) fd) 
(tyo (fix (product outer_var_off_var 1000)) fd) 


(tyo (fix (quotient (product inner_var_off_var 1000) 65536)) £d) 
(tyo (fix (quotient 
[diff (product inner_var_off_var 1000) (product 
(fix (quotient (product inner_var_off_var 1000) 
65536)) 65536)] 256)) £d) 
(tyo (fix (product inner_var_off_var 1000)) fd) 
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s RKAKKKKKKKKKKKKKAKKKK KK KK RKKKK KA RKKKRKKKKA KKK RKKRRKARKAKKARKRKRK 


; READ-DATA * 
; reads data from a file x 
; Inputs: * 
> fd: filename to read data in from * 


/ 
o RARKKKKKRRKKKKKRKRKKKAKARKK KKK KKK KKK RRRKKARARKKKAAAKRKRARARK KKK 


(defun read-data (fd) 
(setq outer_on_road_mean (+ (* 256 (tyi fd)) (tyi fd))) 
(setq inner_on_road_mean (+ (* 256 (tyi fd)) (tyi fd))) 
(setq outer_dev_on_mean (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyi fd))) 1ecom® 
(setq inner_dev_on_mean (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyi £fd))) 100g 


(setq outer_on_road_var (+ (* 256 (tyi fd)) (tyi fd))) 
(setq inner_on_road_var (+ (* 256 (tyi fd)) (tyi fd))) 
(setq outer_dev_on_var (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyi £d))) 10cce 
(setq inner_dev_on_var (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyil £d))) 1ecge 


(setq outer_off_road_mean (+ (* 256 (tyi fd)) (tyi fd))) 
(setq inner_off_road_mean (+ (* 256 (tyi fd)) (tyi fd))) 
(setq outer_dev_off_mean (quotient (float (+ (* 65536 (tyi fd 
))(* 256 (tyl £d)) (tyi td erccn 
(setq inner_dev_off_mean (quotient (float (+ (* 65536 (tyi fd 
)) (* 256 (tyl £d)) (tyr dye ores 


(setq outer_off_road_var (+ (* 256 (tyi fd)) (tyi fd))) 
(setq inner_off_road_var (+ (* 256 (tyi fd)) (tyi fd))) 
(setq outer_dev_off_var (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyi £d))) 1000 
(setq inner_dev_off_var (quotient (float (+ (* 65536 (tyi fd ))(* 
256 (tyi fd )) (tyi £d))) 100g 


(setq outer_mean_threshold (quotient (float(+(* 65536 (tyi fd))(* 
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(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


(setq 


Zoomceyi Lay) Cty L £d)) )e2000))) 
inner_mean_threshold (quotient (float(+(* 65536 (tyi f£d))(* 
256 (tyi fd )) (tyi fd))) 1000)) 
outer_var_threshold (quotient (float(+(* 65536 (tyi fd))(* 
Zoom ita) (tyi fd))) 1000) ) 
inner_var_threshold (quotient (float(+(* 65536 (tyi fd))(* 
Comer) ) (tyi f£d))) 1000) ) 


var_prod_on_road (quotient (float(+(* 65536 (tyi fd))(* 256 
(tyi fd )) tyi fd))) 1000)) 
var_prod_off_road (quotient (float(+(* 65536 (tyi fd))(* 
Z25G Ey ie cadeyecyl £d))) 1000)) 
on_road_prob (quotient (float(+(* 65536 (tyi fd))(* 256 
Gaya fd )) <tyi fd))) 1000)) 
off_road_prob (quotient (float(+(* 65536 (tyi fd))(* 256 
(Gy ec ctv aera) 1000) ) 


inner_on_road_mean (quotient (float(+(* 65536 (tyi fd))(* 
256 (tyi fad )) (tyi fd))) 1000)) 
inner_var_on_mean (float(+(* 65536 (tyi fd))(* 256 (tyi 
fom) tyatd))) ) 
outer_var_on_var (float(+(* 65536 (tyi fd))(* 256 (tyi fd 
)) (tyi fd))) ) 
inner_var_on_var (quotient (float(+(* 65536 (tyi fd))(* 256 
(cya f£diypmtyi £d))) 1000)) 


inner_off_road_mean (quotient (float(+(* 65536 (tyi fd))(* 
ZooeGiy teem) ecky 1 £a))) 1000) ) 
inner_var_off_mean (float(+(* 65536 (tyi fd))(* 256 (tyi 
fd )) (tyl fd))) ) 
outer_var_off_var (quotient (float(+(* 65536 (tyi fd))(* 
25emer werd) je tyierd))) L000) ) 
inner_var_off_var (quotient (float(+(* 65536 (tyi fd))(* 
ZoemCry mec) Cyte£a))) 1000) ) 
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s KAKKKKKKKKKKKAAKKKKKAKAAKKKKAKAKKKKKAKKKRAKKAKKKKAKKAKKRKKAAKKKKK KKK 


; SHOW-ROAD 

; display road already found 

tne ies 

; filename; file where road is stored 

; offset : x screen coordinate when image starts 
» Functions: 


; read-data 


+ % 2 Ht ete Fe He 


; colorcode_l 
4g KAAAKKAAKAKAKKAKKAKKKRKK KAKA KKK KKK KKK KAKA KKK KKK KKK KARA KKK 


(defun show-road (filename offset) 
(setq fd2 (infile filename) ) 
(read-data fd2) 
(setq size (tyi fd2)) 
(colorcode_1l size) 
(do () 
((= (tyipeek fd2) -1)) 
(setq color (tyi fd2)) 
(setg x (+ (+ (* (tyi £02) 256)(tyi £dZ) Sottcer 
(setq y (+ (* (tyi fd2) 256)(tyi fd2)) ) 
(cond ( (= color 1) (writesquare_l x y size 0 255 0)) 


( (= color 2) (writesquare_1l x y size 0 0 255)) 
( (= color 3) (writesquare_l x y size 255 255 0)) 
( (= color 4) (writesquare_l x y size 0 255 255)) 
( (= color 5) (writesquare_l x y size 255 0 255)) 
( (= color 6) (writesquare_l x y size 150 100 100)) 
( (= color 7) (writesquare_l x y size 255 0 0)) 
) 
) 
(close ade} 
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g KAKKKKAAAAKKKKAARAAKAKAKAKAKKKKAKRKKK AKA AA RK AAA AAR AAAK KAKA AKK A 


; PRINTIMSG 


; prints all values on screen 


x 


* 


g RAKAKAAAAKAKARKAAKKAKKKAKAAAK AK KAKKKRKKRARRAKA RRA RRA AKKAKKR AKA A 


(defun 
(msg "' 
(msg "' 
(msg " 
(msg " 


(msg " 
(msg "' 
(msg "' 
(msg " 


(msg " 
(msg "' 
(msg " 
(msg " 


(msg "' 
(msg " 
(msg " 
(msg "' 


(msg " 
(msg " 
(msg " 
(msg " 


(msg " 
(msg " 
(msg " 
(msg " 


printmsg () 

outer_on_road_mean 
inner_on_road_mean 
outer_dev_on_mean 


inner_dev_on_mean 


outer_on_road_var 
inner_on_road_var 
outer_dev_on_var 


inner_dev_on_var 


outer_off_road_mean 
inner_off_road_mean 
outer_dev_off_mean 


inner_dev_off_mean 


outer_off _road_var 
inner_off_road_var 
outer_dev_off_var 


inner_dev_off_var 


outer_mean_threshold 


inner_mean_threshold 


outer_var_threshold 


inner_var_threshold 


var_prod_on_road 
var_prod_off_road 
on_road_prob 


off_road_prob 


outer_on_road_mean N) 
inner_on_road_mean N) 
outer_dev_on_mean N) 


inner_dev_on_mean N) 


outer_on_road_var N) 
inner_on_road_var N) 
outer_dev_on_var WN) 


inner_dev_on_var  N) 


outer_off_road_mean N) 
inner_off_road_mean N) 
outer_dev_off_mean N) 


inner_dev_off_mean N) 


outer_off_road_var N) 
inner_off_road_var N) 
outer_dev_off_var  N) 


inner_dev_off_var N) 


outer_mean_threshold N) 
inner_mean_threshold N) 
outer_var_threshold N) 


inner_var_threshold N) 


var_prod_on_road N) 
var_prod_off_road N) 
on_road_prob N) 
off_road_prob N) 
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(msg 
(msg 
(msg 
(msg 
(msg 
(msg 
(msg 
(msg 


inner_on_road_mean 
inner_var_on_mean 
outer_var_on_var 
inner_var_on_var 
inner_off_road_mean 
inner_var_off_mean 
outer_var_off_var 


inner_var_off_var 


inner_on_road_mean 
inner_var_on_mean 
outer_var_on_var 
inner_var_on_var 
inner_off_road_mean 
inner_var_off_mean 
outer_var_off_var 


inner_var_off_var 
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N) 
N) 
N) 


N) 
N) 
N) 


[ RRRKKKRRKARARKRARKRER AKER RARAR AKA KARKAAEKAER KARR RERKKERAERKKRR / 
/* this is an IRIS-2400 program */ 
/* this is file digit.c 
Its purpose is to read throught a digitizer file 
and produce an image on the IRIS... 
a / 
#include "gl.h" 
#include "device.h" 
display (filename) 


char filename[{26]; 


int fd; /* file descriptor */ 
static RGBvalue bufferl[4096]; 
long i; /* temp loop counter */ 
short lines, cols, type; 

short mask = 255; 

char comment[80]; 

emer bytl([2], byt2[2], byt3[2]; 
/* init the graphics system */ 
pant (); 

/* put into rgb mode */ 
RGBmode () ; 


/* configure the iris...*/ 


gconfig(); 
/* set the color to cyan */ 
RGBcolor(0,255,255); 


/* draw a cyan background */ 


Bectti(0,0,1023,768); 


/* open the file... */ 
fd = open(filename,0); 


if(fd < 0) 


i 


printf("problem in C file indictes not open! n"); 
exit(1); 


/* read the header information off the file */ 
read(fd,byt1,2); 

read(fd,byt2,2); 

read(fd,byt3,2); 

read(fd,comment,80); 


type = bytl[{1] << 8; 

type = type (mask & byti[0]); 
lines = byt2[1] << 8; 

lines = lines (mask & byt2[0]); 
byes) |) << 6; 

cols (mask & byt3[0]); 


cols 


cols 


/* read through the rows of the file... */ 


for(i=0; i < lines; i=i+l) 


/* read a row from the file */ 


read(fd,bufferl,cols); 


/* draw that line */ 
cmov2i(0,(766-1)); 
writeRGB(cols,bufferl ,bufferl ,bufferl); 


close(fd); /* close the file */ 
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[ RAKKRARKKAKKAKAAKKARKKAAARKKRAAKAARAKKAAKAKKAAKAKAKERAKKAKERAKKAKKAKRK / 


/* windowl This is an IRIS 2400 program 
this sets up the cursor and allows the mouse button 
to be used to pick off points, it returens the lower 
left conor of the cursor. 

/* 

#include "gl.h" © 

#include "device.h" 


#include <stdio.h> 


Windowl(size, xmax, ymax, xcursorpos, ycursorpos, type) 
int size, xmax, ymax, type; 


short int xcursorpos[], ycursorpos[]; 


static unsigned short cursordefl1[16]= 
Oxf£800,0x£800,0xd800,0xf£800 , 

Ox£800 ,0x0000 ,0x0000 ,0x0000, 

0x0000 ,0x0000 ,0x0000,0x0000, 

0x0000 ,0x0000,0x0000,0x0000_ ; 

static unsigned short cursordef2[16]= 
ieeecoo), OxfcO0,Oxcc00,0xcc00 , 
Oxfc00,0xfc00 ,0x0000,0x0000, 

0x0000 ,0x0000,0x0000,0x0000, 
0x0000,0x0000,0x0000,0x0000 ; 

static unsigned short cursordef3[16]= 
Oxfe00,O0xfe00 ,0xc600,0xc600 , 

Oxc600 ,Oxfe00 ,Oxfe00 ,0x0000, 

0x0000 ,Ox0000 ,0x0000,0x0000, 

0x0000 ,0x0000,0x0000,0x0000 =; 

static unsigned short cursordef4[16]= 
Oxff00,Oxff£00,0xc300,0xc300 , 
Oxc300,0xc300,Oxff00,Oxff00, 

0x0000 ,0x0000 ,0x0000 ,0x0000, 
0x0000,0x0000,0x0000,0x0000_ ; 

static unsigned short cursordef5[16]= 
Oxf£80,0xff£80,0xcl80,0xcl180 , 


1) 


0xc180 ,0xc180,0xc180 ,Oxff£80, 
Oxf£80,0x0000 ,0x0000,0x0000, 

0x0000 ,0x0000,0x0000,0x0000_ ; 

static unsigned short cursordef6[16]= 
Oxffc0 ,Oxffc0,O0xc0c0,O0xcOcO , 

OxcO0cO0 ,Oxc0cO ,Oxc0c0,0xcO0cO, 
Oxffc0,0xffc0 ,0x0000,0x0000, 

0x0000 ,0x0000,0x0000,0x0000_ ; 

static unsigned short cursordef7[16]= 
Oxffe0,Oxffe0,0xc060,0xc060 , 

Oxc060 ,0xc060 ,0xc060 ,0xc060, 

Oxc060 ,O0xffe0,O0xffe0,0x0000, 
0x0000,0x0000,0x0000,0x0000_ ; 

static unsigned short cursordef8[16]J= 
Oxfff0,0xff£0,0xc030,0xc030 , 

0xc030 ,0xc030 ,0xc030,0xc030, 
0xc030,0xc030 ,Oxfff0,Oxfff0, 
0x0000,0x0000,0x0000,0x0000_ =; 

static unsigned short cursordef9[16]= 
Ox£fES ,OxFErS, OxcOls Oxcold - 
0xc018,0xc018,0xc018,0xc018, 
0xc018,0x%c018,0xc018 ,Oxfff8, 
Oxf££8,0x0000,0x0000,0x0000_ =; 

static unsigned short cursordef10[16]= 
Oxfffc,Oxfffc,Oxc00c,0xc00c , 

OxcO00c ,Oxc00c,0xc00c,0xc00c, 

OxcO00c ,0xc00c,0xc00c,0xc00c, 
Oxfffc,Oxfffc,0x0000,0x0000- ; 

static unsigned short cursordefl1[16]= 
Oxfffe ,Oxfffe,0xc006,0xc006 , 

O0xc006 ,0xc006 ,0xc006 ,0xc006, 

O0xc006 ,0xc006 ,0xc006,0xc006, 
Oxc006,0xfffe,Oxfffe,0x0000_ ; 

static unsigned short cursordef12[16]= 
Oxffff,Oxffff,O0xc003,0xc003 , 
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0xc003 ,0xc003 ,0xc003,0xc003, 
0xc003,0xc003 ,0xc003,0xc003, 
0xc003,0xc003,0xffff,Oxfffrft ; 


/* value returns form the event queue */ 
short value,valuex,valuey; 


/* device name returned from the event queue */ 


Device temp; 


RGBmode() ; 
RGBwritemask(Oxffff ,Oxffff,OxfffFf) ; 


setvaluator (MOUSEX,xmax,0,xmax) ; 
setvaluator (MOUSEY, ymax, 766-ymax,766) ; 


e 
g 


defcursor(l , cursordefl ) 
defcursor(2 , cursordef2 ); 
defcursor(3 , cursordef3 ); 
defcursor(4 , cursordef4 ); 
defcursor(5 , cursordef5 ); 
defcursor(6 , cursordef6 ); 
defcursor(7 , cursordef7 ) 

) 


defcursor(8 , cursordef8 ); 
defcursor(9 , cursordef9 ); 
defcursor(10, cursordefl10) ; 
defcursor(ll, cursordefl1l); 


defcursor(12, cursordef12); 


/* make button gits on the mouse go to the event queue */ 


qdevice(MOUSE1) ; 


/* if button 1 is hit, report mousex and y */ 
tie (MOUSE1 ,MOUSEX ,MOUSEY) ; 


attachcursor (MOUSEX ,MOUSEY) ; 
curson(); 
curorigin(size,2,2); 


RGBcursor(size,255,0,0,0xffff,Oxffff ,OxfffFf) ; 


uo 


RGBcolor (0,255,255); 
rectf£i(550,400, 1023,768); 
RGBcolor(0,0,0); 
cmov2i(550,500); 
while (getbutton(MOUSE1) != 0); 
if(type == 1) 

charstr("choose onroad position"); 


else charstr(" choose offroad position"); 
greset(); 
while(qtest() == 0); 


while(qtest() != 0 ) 


/* 
printf(" got inside the while nn"); 
a 
if(qread(&value) == MOUSE1) 


cursoff(); 


temp = qread(&valuex) ; 
temp = gqread(&valuey) ; 
RGBcolor (0,255,255); 
rectf£i(550,400, 1023,768); 


xcursorpos[0] = valuex ; 

ycursorpos[0] = valuey ; 
/* clean-up */ 

unqdevice(MOUSE1); 

greset(); 

return(type) ; 
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[ RRKRKKKKKAKKAKAKKKAAKAKKAKARAKKRAKKKAKAKKEKRAKK KARR REKKAAKEREAEK / 


/* Readsquare 
reads a square from the screen withe lover left screen 
coordinates of x,y . A square of size size, and returns 


the red bit plane value in the array rbyte 
/ 


mane tude “gl .h" 
#include "device.h'' 


#Hinclude <stdio.h> 


readsquare(x, y, size, rbyte) 
int x, y, size; 


short int rbyte[]; 


int oi: /% x,Y x) 
RGBvalue r[4096], g[4096], b[4096]; 
RGBmode(); 


for (j = 0; j < size; j+t) 


cmov2i(x, ytj); 
readRGB(size, r, g, b); 
for (i = 0; i < size; i++) 


rbyte[j * size + i] = r[i]; 
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[ KRRKRKRKKRKRKRKAKKKKARKARK AKA KKK RKRKK KKK RKKAKKRKRAKKAKKKKK | 


/* Writesquare 
writes a square of size size, at screen coordinates 


x, y. In the mapcolor found with red, green, and blue 


7, 


writesquare(x, y, size, red, green, blue) 


int x, y, size, red, green, blue; 


RGBmode() ; 
RGBcolor(red, green, blue); 


rectfi(x,y, xtsize, ytsize); 
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[RRRKKKKAKKKKKKARKAKAERKAAKKARKKK KAKA KARA REAR ARAKKAKKRK / 
/* Colorcode 
Prints on the screen the color codes used for the seven 
color scan. 
=/ 
colorcode(size) 


int size; 


RGBmode () ; 


RGBcolor(0, 0, 0); 
rectfi(5, 90, 400, 240); 


RGBcolor(0, 255, 0); 
rectfi(20, 220, 20+size-1, 220+size-1); 
cmov2i(40, 220); 


charstr("mean'"') ; 


RGBcolor(0, 0, 255); 
rectfi(20, 200, 20+size-1, 200+size-1); 
cmov2i(40, 200); 


charstr("inner variance"); 


RGBcolor(255, 255, 0); 
rectfi(20, 180, 20+size-1, 180+size-1); 
cmov2i(40, 180); 


charstr("outer variance"); 


RGBcolor(0, 255, 255); 
rectfi(20, 160, 20+size-1, 160+size-1); 
cmov2i(40, 160); 


charstr("mean and inner variance"); 


RGecolor(255, 0, 255); 
rectf1(20, 140, 20+size-1, 140+size-1); 
cmov2i(40, 140); 


charstr("both variance"); 


RGBcolor(150, 100, 100); 
rectfi(20, 120, 20+size-1, 120+size-1); 
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cmov2i(40, 120); 

charstr("mean and outer variance"); 
RGBcolor(255, 0, 0); 

rectfi(20, 100, 20+size-1, 100+size-1); 
cmov2i(40, 100); 


charstr("mean and both variance’); 
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10. 


Pl. 
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